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ABSTRACT 

We summarize the first results from the Gould belt survey, obtained toward the Aquila Rift and Polaris Flare regions during the 'science demon- 
stration phase' of Herschel. Our 70-500 /jm images taken in parallel mode with the SPIRE and PACS cameras reveal a wealth of filamentary 
structure, as well as numerous dense cores embedded in the filaments. Between ~ 350 and 500 prestellar cores and ~ 45-60 Class protostars 
can be identified in the Aquila field, while ~ 300 unbound starless cores and no protostars are observed in the Polaris field. The prestellar core 
mass function (CMF) derived for the Aquila region bears a strong resemblance to the stellar initial mass function (IMF), already confirming the 
close connection between the CMF and the IMF with much better statistics than earlier studies. Comparing and contrasting our Herschel results in 
Aquila and Polaris, we propose an observationally-driven scenario for core formation according to which complex networks of long, thin filaments 
form first within molecular clouds, and then the densest filaments fragment into a number of prestellar cores via gravitational instability. 

Key words, stars: formation - stars: luminosity function, mass function - stars: protostars - ISM: clouds - ISM: structure - submillimeter: ISM 



1. Introduction 

The Herschel Space Observatory (Pilbratt et al. 2010) offers a 
unique opportunity to improve our global understanding of the 
earliest phases of star formation. Here, we present first highlights 
from the Gould belt survey, one of the largest key projects with 
Herschel (cf. Andre & Saraceno 2005), based on extensive far- 
infrared and submillimeter mapping of nearby molecular clouds 
with both the SPIRE (Griffin et al. 2010) and PACS (Poglitsch et 
al. 2010) bolometer cameras. This SPIRE/PACS imaging survey 
will cover the bulk of the nearest {d < 0.5 kpc) cloud complexes 
in the Galaxy, which are mostly located in the Gould belt, a gi- 
ant (~ 700 pc X 1000 pc), flat structure incHned by ~ 20° to the 
Galactic plane (e.g., Guillout 2001). Since the ~ 15" angular 
resolution of Herschel around A ~ 200 //m is adequate for prob- 
ing individual (~ 0.01-0.1 pc) star-forming cores up to ~ 0.5 kpc 
away, the cloud complexes of the Gould belt correspond to the 
volume of Galactic space where Herschel imaging can be best 
used to characterize in detail the earliest stages of star formation. 
The immediate observational objective of the Gould belt sur- 
vey is to obtain complete samples of nearby prestellar cores and 
Class protostars with well characterized luminosities, tempera- 
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lures, and density profiles, as well as robust core mass functions 
and protostar luminosity functions, in a variety of star-forming 
environments. An order of magnitude more cold prestellar cores 
than already identified from the ground are expected to be found 
in the entire survey, which should allow us to derive an accu- 
rate prestellar core mass function (CMF) from the pre-brown- 
dwarf to the intermediate-mass range. Thanks to its high sen- 
sitivity and large spatial dynamic range, this Herschel survey 
can also probe, for the first time, the link between low-density 
cirrus-like structures in the interstellar medium (ISM) and com- 
pact self-gravitating cores. The main scientific goal is to eluci- 
date the physical mechanisms responsible for the formation of 
prestellar cores out of the diffuse ISM, which is crucial for un- 
derstanding the origin of the stellar initial mass function (IMF). 
Our first results, obtained toward the Aquila Rift and Polaris 
clouds, are very promising (e.g., Konyves et al. 2010, Bontemps 
et al. 2010, Men'shchikov et al. 2010, Ward-Thompson et al. 
2010, and Miville-Deschenes et al. 2010). As discussed in 
Sect. 4 of this paper, they suggest that prestellar cores result from 
the gravitational fragmentation of filaments in the cold ISM. 



2. Herschel observations 

The Herschel survey was designed to cover the densest portions 
of the Gould belt with SPIRE at 250-500 yum and PACS at 100- 
160 fim. The observational goal is to make a complete, homo- 
geneous mapping of the Ay > 3 regions with SPIRE and the 
Ay > 6 regions with PACS, and representative areas at Ay ~ 1-3 
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levels with both instraments. The survey sensitivity (better than 
Ay ~ 1 at the 5cr level) will allow us to probe the structure of 
nearby molecular clouds down to the level of the interface with 
their atomic gas envelopes. The 15 target clouds span a range of 
physical conditions, from active, cluster-forming complexes to 
quiescent regions with lower or no star formation activitjQ. 
Our Herschel mapping consists of two steps: 

(1) A wide-field SPIRE/PACS survey of a total surface area 
~ 160 deg^ using the so-called "parallel mode" with a scan- 
ning speed of 60"/s and the PACS 70 /im and 160 yum bands. 
In this first step, the main goal is to acquire adequate SPIRE 
250-500 yum data for all of the target regions. The PACS data 
acquired simultaneously yield very useful information at 70 /im 
and 160 /^m through most of the SPIRE survey, but do not have 
optimal (diffraction-limited) angular resolution. We selected the 
70 jim filter for the blue band of PACS to obtain a good diag- 
nostic of the presence of embedded protostars (cf. Dunham et al. 
2008) throughout the mapped regions. 

(2) A dedicated PACS-only survey of a total surface area 
~ 65 deg^, observing the 100 jivcl and 160 //m bands with a scan- 
ning speed of 20"/s. This second step is supposed to provide data 
of optimal sensitivity and resolution at 100 /im and 160 /im. 

The Aquila Rift and Polaris flare regions were observed in 
parallel mode with both SPIRE and PACS (step 1) during the 
science demonstration phase of Herschel. The coiTesponding 
observations are described in detail by Konyves et al. (2010), 
Bontemps et al. (2010), Men'shchikov et al. (2010), and Ward- 
Thompson et al. (2010) in this special A&A issue. The Polaris 
flare field is a high-latitude translucent cloud with little to no star 
formation ?Ad ~ 150 pc (e.g., Heithausen et al. 2002), and is ex- 
pected to have the lowest level of background cloud emission 
and cirrus confusion noise in the entire Gould belt survey. At 
the other extreme, the Aquila field is a very active star-forming 
complex at c/ ~ 260 pc (e.g., Gutermuth et al. 2008), and is 
expected to have the highest level of background cloud emis- 
sion and cirrus confusion noise in the whole survey. Sensitive 
Galactic far-IR imaging surveys such as the one discussed here 
are not limited by instrumental sensitivity but by confusion aris- 
ing from small-scale cirrus/cloud structure (cf. Roy et al. 2010 
and references therein). The levels of background cloud fluc- 
tuations observed in the Aquila and Polaris fields allow us to 
roughly estimate the whole range of ciiTus confusion noise lev- 
els that will affect the Gould belt survey. In the Aquila field, the 
observed rms level of cirrus or "structure" noise (as measured 
on the typical scales of dense cores) ranges from ~ 30 mJy/18"- 
beam to ~ 300 mJy/18"-beam at /I = 250 yum, which is a factor 
of ~ 3-30 above the rms instrumental noise. In contrast, in the 
Polaris field, the rms level of emission fluctuations measured in 
the SPIRE 250jum map is only ~ 10-30 mJy/18"-beam, which is 
very close to the rms sensitivity level expected after two cross- 
scans in parallel mode at 60"/s. Assuming the dust opacity of 
Hildebrand (1983) and the median dust temperatures derived 
from our Herschel maps, the above surface-brightness sensitiv- 
ity levels translate into 5cr column-density detection thresholds 
of A^Hl ^ 3 X 10^° cm"^ in Polaris and A^h2 ^ 10^' cm"^ in 
Aquila (in agreement with the column-density maps shown by 
Ward-Thompson et al. 2010 and Konyves et al. 2010). This high 
column-density sensitivity, coupled with the unprecedented sur- 
face brightness and spatial dynamic ranges of our Herschel map- 
ping, allows us to probe, for the first time, the physical connec- 
tion between the structure of the diffuse ISM and the formation 
of prestellar cores (cf. Sect. 4). 
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The corresponding mass sensitivities for typical prestel- 
lar/protostellar cores are ~ 0.3 Mq (85% completeness level) in 
Aquila {d - 260 pc) and ~ 0.01 Mq (85% completeness level) 
in Polaris (d - 150 pc). These completeness numbers were es- 
timated by performing Monte-Carlo simulations as described by 
Konyves et al. (2010) in the case of Aquila. The initial results 
obtained in Aquila and Polaris therefore confirm that the com- 
pleteness level of our Herschel census for prestellar cores will 
reach the pre -brown dwarf mass regime (Mcoie < 0.08 Mq) in the 
nearest molecular clouds (d ~ 150 pc) of the Gould belt. 



3. Main results and analysis 

The Herschel images of the Aquila Rift and Polaris flare re- 
gions exhibit extensive filamentary structure, as well as numer- 
ous dense cores situated along these filaments (see Fig. 1, on- 
line Fig. 3, and Men'shchikov et al. 2010). A total of 541 star- 
less cores (~ 0.01-0.1 pc in size) can be identified in the whole 
(> 3.3° X 3.3°) Aquila field, most (> 60%) of which appear 
to be self-gravitating and prestellar in nature. The latter is in- 
ferred from: (1) a comparison of the core masses derived from 
the SPIRE/PACS spectral energy distributions with local values 
of the Jeans or Bonnor-Ebert (BE) mass also estimated from 
Herschel data (see Konyves et al. 2010); (2) the positions of the 
cores in a mass versus size diagram (cf. online Fig. 4), which are 
close to the loci expected for critical isothermal BE spheres at 
gas temperatures ~ 7-20 K; and (3) the mean column densities 
of the cores, which exceed the background column densities by 
a median factor ~ 1 .5 as expected for critically self-gravitating 
BE spheres. The shape of the core mass function (CMF) de- 
rived for this sample of Aquila cores closely resembles the IMF 
(Fig. 2-left - see Konyves et al. 2010 for the CMF of the Aquila 
central region and Appendix A for details on the derivation of 
core masses). In contrast, the 302 starless cores identified in the 
ciiTus-like Polaris cloud with the same clump-finding algorithm 
igetsources - see Men'shchikov et al. 2010) appear to be mostly 
unbound (cf. online Fig. 4) and their mass distribution does not 
follow the IMF, with a peak at an order of magnitude smaller 
mass (Fig. 2-right). Only 5 of the Polaris cores are reasonably 
close to being gravitationally bound and thus possibly prestellar 
in nature (see Ward-Thompson et al. 2010). Between 45 and 60 
Class protostars are revealed by Herschel in the Aquila field 
(Bontemps et al. 2010), while not a single protostar is detected 
in the Polaris region. 

The Aquila filaments harboring embedded protostars and/or 
large concentrations of prestellar cores are all characterized by 
higher column densities (A^H, ~ 10^^ cm"^), suggesting that they 
are gravitationally unstable (see Fig. 1-left). In contrast, both the 
Polaris filaments and the quiescent, non-star-forming filaments 
observed in Aquila have much lower column densities (up to a 
few 10^' cm"-), suggesting they are gravitationally stable (see 
Fig. 1-right and Sect. 4 below). 



4. Discussion and conclusions 

The Herschel early results obtained toward the Aquila star- 
forming complex confirm that the prestellar CMF resembles the 
stellar IMF, using data with already a factor of ~ 2 to 9 bet- 
ter counting statistics than earlier (sub)-millimeter continuum or 
near-IR extinction surveys (e.g., Motte et al. 1998, Johnstone 
et al. 2000, Enoch et al. 2006, Stanke et al. 2006, Nutter & 
Ward-Thompson 2007, Alves et al. 2007). The close resem- 
blance between the Aquila CMF and the IMF in both shape and 
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Fig. 1. Column density maps of two subfields in Aquila (left) and Polaris (right) derived from our SPIRE/PACS data. The contrast 
of the filaments with respect to the non-filamentary background has been enhanced using a curvelet transform as described in 
Appendix A. Given the typical width ~ 10000 AU of the filaments, these column density maps are equivalent to maps of the mass 
per unit length along the filaments. The color scale shown on the right of each panel is given in approximate units of the critical line 
mass of Inutsuka & Miyama (1997) as discussed in Sect. 4. The areas where the filaments have a mass per unit length larger than 
half the critical value and are thus likely gravitationally unstable have been highlighted in white. The maximum line mass observed 
in the Polaris region is only ~ 0.45 x the critical value, suggesting that the Polaris filaments are stable and unable to form stars at 
the present time. The candidate Class protostars and bound prestellar cores identified in Aquila by Bontemps et al. (2010) and 
Konyves et al. (2010) are shown as green stars and blue triangles, respectively. Note the good correspondence between the spatial 
distribution of the bound cores/protostars and the regions where the filaments are unstable to gravitational collapse. 



mass scale suggests that, as a class, the self-gravitating prestel- 
lar cores identified in far-IR/submilimeter continuum imaging 
studies such as the present Herschel survey may form stars on a 
one-to-one basis, with a fixed and relatively high local efficiency, 
i.e., fcoi-e = Mi,/Mcore ~ 20 - 40% in Aquila. This is consistent 
with theoretical models according to which the stellar IMF is in 
large part determined by pre-collapse cloud fragmentation, prior 
to the protostellar accretion phase (cf. Larson 1985; Padoan & 
Nordlund 2002; Hennebelle & Chabrier 2008). There are sev- 
eral caveats to this simple picture (cf. discussion in Andre et al. 
2009), and detailed analysis of the data from the whole Gould 
belt survey will be required to fully characterize the CMF-IMF 
relationship and, e.g., investigate possible variations in the ef- 
ficiency fcore with environment. It is nevertheless already clear 
that one of the keys to the problem of the origin of the IMF lies 
in a good understanding of the formation process of prestellar 
cores, even if additional processes, such as rotational subfrag- 
mentation of prestellar cores into binary /multiple systems (e.g.. 
Bate et al. 2003), probably also play an important role. 

Our Herschel initial results also provide key insight into the 
core formation issue. They support an emerging picture (see also 
Myers 2009) according to which complex networks of long, thin 
filaments form first within molecular clouds, possibly as a result 
of interstellar MHD turbulence, and then the densest filaments 
fragment into a number of prestellar cores via gravitational in- 
stability. That the formation of filaments in the diffuse ISM rep- 
resents the first step toward core/star formation is suggested by 
the filaments already being omnipresent in a diffuse, non-star- 
forming cloud such as Polaris (cf. Fig. 1-right, Men'shchikov et 



al. 2010, and Miville-Deschenes et al. 2010). The second step 
appears to be the gravitational fragmentation of a subset of the 
filaments into self-gravitating cores. This observationally-driven 
scenario can be placed on a stronger footing by comparing our 
Herschel results to existing models of filamentary cloud frag- 
mentation. Inutsuka & Miyama (1992, 1997) showed that an 
unmagnetized isothermal filament is unstable to axisymmetric 
perturbations if the line mass or mass per unit length, Mune, of 
the filament is larger than the critical value required for equi- 
librium, M."""''^, = Ic'JG, where c. is the isothermal sound 

' line.cnt .? ' ' ^ 

speed. Remarkably, the critical line mass only depends on gas 
temperature (Ostriker 1964) and is modified by only a factor 
of order unity for filaments with realistic levels of magnetiza- 
tion (Fiege & Pudritz 2000). Figure 1 shows maps of the mass 
per unit length, expressed in approximate units of the critical 
line mass, for the filaments detected by Herschel in Aquila and 
Polaris. These maps were constructed from the column density 
maps derived from our SPIRE/PACS images (see Appendix A) 
by multiplying the local column density of the filaments by their 
measured typical width (FWHM ~ 14000 AU in Aquila and 
~ 9000 AU in Polaris, before deconvolution - cf. Men'shchikov 
et al. 2010). In Fig. 1, the critical line mass coiTesponding to a 
gas temperature of 10 K (A^ii,-,e t™™'*^ ~ 15 Mo/pc) was adopted 
throughout the fields, but nearly identical results are obtained 
if the dust temperature maps derived from the Herschel im- 
ages (cf. Bontemps et al. 2010) are used instead. The results (cf. 
Fig. 1) show that most (> 60%) of the bound prestellar cores 
identified in Aquila are concentrated in supercritical filaments 

withMH„e>M;;:^, 



Furthermore, virtually all supercritical fil- 
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Fig. 2. Core mass functions (blue histograms with error bars) derived from our SPIRE/PACS observations of the Aquila (left) and 
Polaris (right) regions, which reveal of total of 541 candidate prestellar cores and 302 starless cores, respectively. A lognormal fit 
(red curve) and a power-law fit (black solid line) to the high-mass end of the Aquila CMF are superimposed in the left panel. The 
power-law fit has a slope of -1.5 + 0.2 (compared to a Salpeter slope of -1.35 in this dA^/dlogM format), while the lognormal fit 
peaks at ~ 0.6 M^ and has a standard deviation of ~ 0.43 in logioM. The IMF of single stars (corrected for binaries - e.g., Kroupa 
2001), the IMF of multiple systems (e.g., Chabrier 2005), and the typical mass spectrum of CO clumps (e.g., Kramer et al. 1998) are 
also shown for comparison. Note the remarkable similarity between the Aquila CMF and the stellar IMF, suggesting a ~ one-to-one 
correspondence between core mass and star/system mass with M*sys - eA^coie and e x; 0.4 in Aquila. 



aments harbor prestellar cores and/or embedded protostars, in 
agreement with the view that they are collapsing and actively 
forming stars at the present time. In particular, this is the case 
for the Aquila main filament, which has Mii„e > 5 x A^i"™^fit and 
is associated with a very rich protocluster (Serpens South - see 
Gutermuth et al. 2008, Bontemps et al. 2010). In contrast, the 
subcritical filaments with Miine < -^^iln^rit are generally devoid 
of prestellar cores and protostars, which is consistent with the 
view that they are gravitationally stable, hence neither collaps- 
ing nor forming stars. All of the Polaris filaments appear to be 
in the subcritical regime, the maximum observed value of the 
stability parameter being A/iine/A^i"™"rit ~ 0-45, and it is unclear 
whether they will evolve into the unstable regime or not. 

It is noteworthy that the critical line mass approximately cor- 
responds to a critical column density A^H,,crit ~ 10^^ cm"^, i.e., 
to a critical visual extinction Av,crit ~ 10. Our Herschel findings 
thus provide an explanation of the visual extinction threshold 
for core formation found by earlier ground-based studies (e.g., 
Onishi et al. 1998, Johnstone et al. 2004). Prestellar cores are 
only observed above a threshold Ay.ci-it because they form out of 
a filamentary background and only the supercritical, gravitation- 
ally unstable filaments are able to fragment into bound cores. 

Confirming and refining this scenario for core formation will 
require the results of the entire Herschel survey, as well as 
follow-up (sub-)millimeter dust polarimetry and molecular line 
observations to, e.g., clarify the roles of magnetic fields, turbu- 
lence, and gravity in forming the filaments. Our initial results 
are nevertheless extremely encouraging. Extrapolating from the 
^ 500 prestellar cores and t 45 Class protostars identified 
in the ~ 19 deg^ covered by the Aquila and Polaris fields, we 
expect that the -160 deg^ Gould belt survey will reveal a to- 
tal of about 4500 prestellar cores, including a large number of 
candidate pre-brown dwarfs in the nearest {d ~ 150 pc) clouds, 
and more than 350 Class protostars. This will provide a unique 
database, including the southern hemisphere, for follow-up high- 



resolution molecular line/dust continuum studies of the physics 
of individual cores and protostars with ALMA. 
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Appendix A: Derivation of core/filament properties 
and effects of distance uncertainties 

As described in more detail in a companion paper by Konyves 
et al. (2010) on Aquila, the masses of the cores identi- 
fied in the Herschel images with the getsources algorithm 
(see Men'shchikov et al. 2010) were derived by fitting grey- 
body functions to the spectral energy distributions (SEDs) 
constructed from the integrated flux densities measured with 
SPIRE/PACS for each core. We assumed the dust opacity law 
Ky - 0.1 (v/1000 GHz)^ cmVg, where v denotes frequency 
and Ky is the dust opacity per unit (gas + dust) mass column 
density. This dust opacity law, which is very similar to that 
advocated by Hildebrand (1983), is consistent with the value 
Ki.-imm - 0.005 cm^/g adopted for starless cores in numerous 
earlier studies (e.g., Motte et al. 1998). Ignoring any systematic 
distance effect (see below), the core mass uncertainties are dom- 
inated by the uncertainty in Ky, typically a factor of ~ 2. Cores 
were classified as either protostellar or starless based on the pres- 
ence or absence of significant PACS emission at 70 /im, respec- 
tively (cf. Bontemps et al. 2010 and Dunham et al. 2008). In the 
Polaris field, the cirrus noise level is so low (cf. Fig. 4-right) 
that we cannot exclude that a fraction of the 302 candidate star- 
less cores extracted with getsources correspond to background 
galaxies. 

A column density map was derived for each region from 
the Herschel images smoothed to the SPIRE 500 //m resolution 
(36.9" FWHM) using a similar SED fitting procedure on a pixel 
by pixel basis (see, e.g.. Figs. 1 and 6 of Konyves et al. 2010 
for Aquila). To obtain the maps of the filamentary background 
shown in Fig. 1 of this paper, we then performed a "morphologi- 
cal component analysis" decomposition (e.g., Starck et al. 2003) 
of the original column density maps on curvelets and wavelets. 
The curvelet component images shown in Fig. 1 provide a good 
measurement of the column density distribution of the filamen- 
tary background after subtraction of the compact sources/cores 
since the latter are contained in the wavelet component. We es- 
timate that these column density maps are accurate to within a 
factor of ~ 2. The scaling in terms of the mass per unit length 
along the filaments is more uncertain, however, as it depends on 
distance (see below) and would in principle require a detailed 
analysis of the radial profiles of the filaments, which is beyond 
the scope of the present letter. Here, we simply assumed that 
the filaments had a Gaussian radial column density profile and 

multiplied the surface density maps by Jj^ x W ~ 1.06 W, 

where W is the typical FWHM width of the filaments. We as- 
sumed a mean molecular weight of jj. - 2.33. At this stage, the 
correspondence between the critical line mass of the filaments, 
^unmdg ^^j j.j^^ visual extinction threshold, Avrrit ~ 10 (see 

lme,cnt' ' v,t.iiL v 

Sect. 4), is thus accurate to at best a factor of ^ 2. 

There is some ambiguity concerning the distance to the 
Aquila Rift region. A number of arguments, presented in a com- 
panion paper by Bontemps et al. (2010), suggest that the whole 
region corresponds to a coherent cloud complex at d- - 260 pc 
(see also Gutermuth et al. 2008), which is the default distance 
adopted in the present paper for Aquila. However, other studies 
in the literature (see references in Bontemps et al. 2010) place 
the complex at a larger distance, d+ - 400 pc. It is thus worth 
discussing briefly how our Aquila results would be affected if we 
adopted the larger distance estimate, d+, instead of d-. The core 
mass estimates, which scale as Syd^ l[By{TA)Ky] where Sy is in- 
tegrated flux density and By{Ti) is the Planck function, would 
systematically increase by a factor of 2.4. This would shift the 



CMF shown in Fig. 2-left to the right and thus lower the effi- 
ciency ecoie from ~ 20-40% to ~ 10-20%. In the mass versus 
size diagram of Fig. 4, the cores would move upward as indi- 
cated in the left panel of the figure, which would increase the 
fraction of candidate bound cores in Aquila from 63% to 81%. 
The column density map of the Aquila filaments shown in Fig. 1 
would be unchanged, but the scaling in terms of the mass per 
unit length along the filaments would change by ~ 50% up- 
ward, since the physical width of the filaments would increase 
by ~ 50%. In other words, the highlighted regions in Fig. 1- 
left, where the mass per unit length of the filaments exceeds half 
the critical value, would slightly expand, increasing the contrast 
with the Polaris filaments and improving the correspondence be- 
tween the spatial distribution of the prestellar cores/protostars 
in Aquila and that of the gravitationally unstable filaments. To 
summarize, our main conclusions do not depend strongly on the 
adopted distance. 
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Fig. 3. Composite 3-color images of the ~ 11 deg^ Aquila field (left) and the ~ 8 deg^ Polaris field (right) produced from our 
PACS/SPIRE parallel-mode images at 70, 160, and 500 ;um. The color coding of the Aquila image is such that Red = SPIRE 500 
yum. Green - PACS 160 yum. Blue = PACS 70 fim. For the Polaris image. Red = combination of the three SPIRE bands. Green = 
PACS 160 |/m. Blue = PACS 70 /im. The Aquila composite image was also the first release of "OSHI", ESA's Online Showcase of 
Herschel Images (cf. http://oshi.esa.int and http://www.esa.int/SPECIALS/Herschel/SEMT0T9K73Gi).html). 
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Fig. 4. Mass vs. size diagrams for the starless cores detected with //eric/ie/-SPIRE/PACS in Aquila (left) and Polaris (right) (blue 
triangles). The masses were derived as explained in Appendix A (see also Konyves et al. 2010) and the sizes were measured at 
250 yum. For reference, the locations of the (sub)mm continuum prestellar cores identified by Motte et al. (1998, 2001) in p Oph and 
NGC2068/2071, respectively, are shown, along with the correlation observed for diffuse CO clumps (shaded band - cf. Elmegreen 
& Falgarone 1996). The two black solid lines mark the loci of critically self-gravitating isothermal Bonnor-Ebert spheres at T = 7 K 
and r = 20 K, respectively. The 341 cores classified as prestellar by Konyves et al. (2010), out of a total of 541 starless cores 
in Aquila, are shown as filled triangles in the left panel. The 302 starless cores of Polaris lie much below the two Bonnor-Ebert 
lines in the right panel, suggesting they all are unbound. The red solid lines are two lines of constant mean column density at 
Nh^ - 10^' cm"^ and Nn^ - 10^^ cm~^, respectively. The 5cr detection threshold at li = 150 pc of existing ground-based (sub)mm 
(e.g., SCUBA) surveys as a function of size is shown by the dashed curve. (The shape of this curve reflects a constant sensitivity to 
column density until source size approaches the beam size.) The 5cr detection thresholds of our SPIRE 250 /im observations, given 
the estimated levels of cirrus noise (cf. Sect. 2) and assumed distances of 260 pc for Aquila and 150 pc for Polaris, are shown by 
the blue curves. The blue arrow in the left panel indicates how the positions of the Herschel cores (and the blue curve marking the 
detection threshold) would move in the diagram if a distance of 400 pc were adopted for Aquila instead of 260 pc (see Appendix A). 



